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Abstract. To study the onset of a stationary dynamo in the presence of inner or outer walls of various 
electromagnetic properties, we propose a simple ID-model in which the flow is replaced by an alpha effect. 
The equation of dispersion of the problem is derived analytically. It is solved numerically for walls of 
different thicknesses and of electric conductivity and magnetic permeability different from those of the 
fluid in motion. We also consider walls in the limit of infinite conductivity or permeability. 
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1 Introduction 

A number of experimental devices have been built in the 
last years aiming at producing dynamo action (for reviews 
see e.g. [JJ and Such a device is generally made of 
a container in which some liquid metal is put into mo- 
tion. In a previous study [3] we considered the influence of 
the electromagnetic properties of the container outer wall 
onto the onset of dynamo action for the Riga (Latvia) and 
Karlsruhe (Germany) experiments. The results depend on 
which type of dynamo instability is obtained. For station- 
ary solutions like in the Karlsruhe experiment [4] [5], the 
reduction of the dynamo instability threshold is mono- 
tonic versus the conductivity and the permeability of the 
outer wall. For oscillatory solutions like in the Riga experi- 
ment [5] [7] , there are additional eddy currents in the outer 
wall. These currents produce an additional dissipation op- 
posed to the reduction of the threshold. In that case, the 
reduction of the dynamo instability threshold versus the 
conductivity and the permeability of the outer wall is not 
monotonic anymore. 

These results are consistent with other studies aiming at 
studying the influence of the thickness of a stagnant outer 
layer conducting fluid (or equivalently of an outer wall 
with the same conductivity as the fluid) on the dynamo 
threshold. Various inner flow geometries have been con- 
sidered leading to either stationary [S][S] or oscillatory 
[TO] [IT] solutions. In Fig. [TJ we give a synthesis of these 
results. For that we plot the threshold reduction rate r 
versus e/R where F = 1 — R m {e/R)/R m {Q), e is the thick- 
ness of the stagnant outer layer, R the radius of the fluid 
container and R m the magnetic Reynolds number defined 



email: Franck.Plunian@hmg.inpg.fr 



by R m — UR/rj where U is a characteristic flow inten- 
sity and r\ the magnetic diffusivity. The stationary solu- 
tions correspond to the full curves (a-d) and are increas- 
ing monotonically versus e/R. The non-stationary solu- 
tions correspond to the dashed curves (e-i) and reach a 
maximum versus e/R (though not obvious for the curves 
(f) and (h) it is actually the case). The influence of the 
electric conductivity of an inner core has also been inves- 
tigated [12]. Again the same conclusions as [3] have been 
found: for stationary (resp. oscillating) solutions, the dy- 
namo threshold decreases monotonically (resp. reaches a 
minimum) when increasing the conductivity of the inner 
core. 

Given the high difficulties to build a dynamo experiment, 
reducing the threshold by changing the electromagnetic 
boundary conditions is of course of high interest for the 
experimenter. The influence of electromagnetic boundary 
conditions onto dynamo action has also been the object 
of different studies relevant to planets, and stars. In the 
case of Earth-like planets for example, the influence of a 
conducting solid inner core onto the dynamo action pro- 
duced by the outer core fluid motion has been studied by 
different ways, using either a prescribed a 2 -effect [12], a 
prescribed a-effect and buoyancy [13] or a direct resolu- 
tion of the full convective dynamo model [14j . The main 
issue of these studies was to identify whether the inner 
core has a stabilizing effect on the reversals of the dipole 
component of the magnetic field. So far there is no defi- 
nite answer to this problem as these different studies lead 
to contradictory results. In the case of Solar-like stars, it 
has been shown [15] how some magnetic features (like the 
PDF of the magnetic field strength) observed at the sur- 
face of the star could give some indications on the relevant 
magnetic boundary conditions of a turbulent convective 
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Fig. 1. Threshold reduction rate r versus e/R for (a) the Von 
Karman sodium experiment [8], the Kumar and Roberts flow 
for (b) R m > and (c) R m < [5], (d) the Karlsruhe experi- 
ment [3], (e) the Perm experiment [11], (f) the Riga experiment 
[3], the sltl flow with e = (g) 0.35, (h) 0.3, (i) 0.25 [TO]. 



dynamo model. 

Magnetic boundary conditions may also be important in 
the case of Fast Breeder Reactors (FBR). Indeed, though 
these industrial installations have not being designed to 
produce dynamo action, they share some common features 
with the Karlsruhe experiment. In the core of a FBR, the 
liquid sodium flows in an array of a large number of par- 
allel straight tubes called assemblies. In each of them the 
flow geometry is again composed of a periodic array of 
single helical vortices. One has shown [16] the existence 
of an a-effect similar to the Karlsruhe experiment in each 
assembly. Though such an a-effect is not sufficient to gen- 
erate a dynamo instability for a core with homogeneous 
electromagnetic properties, the question remains when the 
walls of the assemblies are made, for example, of ferro- 
magnetic steel (relative permeability of the order 10 3 and 
relative conductivity of order 1) [17] and when the array 
of assemblies is surrounded by a belt of ferromagnetic ma- 
terial (as it is the case in the FBR Phenix). In that case 
both inner and outer walls electromagnetic properties may 
be important. Indeed a reduction of the dynamo thresh- 
old leading to some dynamo instability within the core of 
the FBR could imped the right working of the reactor and 
lead for example to an emergency breakdown. 

In the present paper we consider a simple ID-model 
of a stationary dynamo in order to evaluate the relative 
importance of inner and outer wall electromagnetic prop- 
erties onto the onset of dynamo action. For that we con- 
sider two types of boundary conditions: either an outer 



wall and isolating medium outside or periodic inner walls. 
In both cases we not only vary the relative conductivity 
and permeability of the walls but also their thicknesses. 



2 A simple ID-model 

2.1 An anisotropic a-effect 

We use the kinematic approach consisting in solving the 
induction equation for a given motion. This equation reads 



dB 

~dt 



Vx(UxB) + i)V 2 B, V.B = 0, 



(1) 



where, again, rj means the magnetic diffusivity of the fluid, 
B the magnetic field and U the fluid velocity. 
As we are not interested by a flow geometry in particular 
but only by the influence of the boundary conditions onto 
the onset of dynamo action, we assume that the interac- 
tions of the flow with the magnetic field can be represented 
by an anisotropic a-effect like in the Karlsruhe experiment 
fTS] [TO] or in the core of a FBR [IB]. Following the lines 
of mean-field dynamo theory [20] the magnetic field B 
and the fluid velocity U are expressed as sums of mean 
fields, B and U, and fluctuating fields, B' and U'. Here 
the mean is defined by a space average of the original field. 
Assuming U = 0, the mean part of the induction equation 
([I]) writes 



dB V7 C 



V V 2 B, 



V.B = 0, 



(2) 



where £ is a mean electromotive force due to the fluid 
motion given by 

£ = U' x B'. (3) 

We may consider £ as a functional of U' and B. Let us 
accept the assumption usually adopted in the mean-field 
context that £ in a given point in space and time depends 
on B only via the components of B and their first spatial 
derivatives in this point. This is reasonable for sufficiently 
small variations of B in space and time. For a first approx- 
imation, on which we restrict ourselves here, we consider 
no other contributions to £ than that describing the a- 
effect, that is, we ignore all contributions to £ containing 
derivatives of B. In addition, the a-effect is assumed to 
act in the xy-plane only, where x, y and z are the Carte- 
sian coordinates. This corresponds to a flow for example 
independent of z. Then ([2]) writes in the form 



— = -Vx[a(B-(z.B)z)] 



7yV 2 B, V.B = (4) 



where a is a scalar quantity. Such a model has proved 
to be sufficiently realistic for both cases, the Karlsruhe 
experiment [H][IH] an d the FBR core [TBJ. However for 
the latter case, there is an additional mean flow U along 
the z-direction which affects the model. The corresponding 
discussion is postponed to section [5] 
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Fig. 2. Scheme of the 3 regions of the non periodic problem 
(above) and of the 2 regions of the periodic problem (below). 



2.2 Model parameters 

We consider three regions I (=1, 2 or 3) symmetric with 
respect to the plane x = and infinite in the y and z 
directions. They are defined by their respective size along 
x (±x < R, R < ±sg < R+e and ±x > i?+e), conductivity 
(tTi, <72, <73 ), permeability (/xi, /i2, /X3) and a-effect (cvi = 
a, a2 = 0, CV3 = 0) where a is a steady scalar quantity 
which does not depend on 1, y nor z. 
We also consider two types of boundary conditions in the 
x-direction (see Fig. [2|. The problem in which the region 
3 is an insulator (03 = 0) and extends to infinity with 
B (x) — > when x — > 00 is called the non-periodic problem. 
In that case the region 2 corresponds to an outer wall. 
The problem in which the region 3 does not exists and 
B(x = R + e) = B(.x = —R — e) is called the periodic 
problem. In that case the region 2 corresponds to periodic 
inner walls like the walls of the assemblies of a FBR. 



2.3 Reduction of the basic equations 

The solutions of f4| can be represented as series of Fourier 
modes proportional to exp(ijy + ikz). As a does not de- 
pend neither on y nor z, each (J, fc)-mode is independent 
from each other. As a is steady we may expect solutions 
varying like exp(pt) in time, with the real part of p be- 
ing the growth rate of the magnetic field. In the rest of 
the paper we consider only the mode j = 0, for sake of 
simplicity. Then, for a given fc, we may look for B in the 
form 



a and b being functions of x only. Replacing (|5|) in ([4]), we 
find the following equations 



{rjiA - p)a t - a t bi = 
(rjiA -p)bi - aik 2 a t = 



(6) 
(7) 



where r\\ — 1/aifii, A = d x 2 — k 2 and Z = 1,2,3 (1 = 1, 2) 
for the non-periodic (periodic) problem. 
We can show that there exists two sets of independent 
solutions depending on the parity of a\(x) and b\(x). In- 
deed, from ([6]) and ||7) ai is solution of L(a%) — with 
L = (ijiA—p) 2 — a 2 k 2 . This operator L is linear and leaves 
the parity of the function unchanged. Therefore writing a\ 
as the sum of an odd and even functions a\ — a° + a\ we 
find that L(a°) = L(a\) = 0. Consequently, a° and a\ are 
two independent solutions of L. From ([6]) and (O, it is 
easy to show that b\ has the same parity as a±. There- 
fore it is sufficient to solve ^ and ([7]) for each parity, the 
general solution being a linear combination of them. The 
solution at x — is given by b = (0, 0, b z ) if a\(x) is odd 
and b = (0, b y ,b z ) if a\(x) is even. 



3 Method of resolution 
3.1 General solutions 

The general form of the solutions of equations ([6|) and ([7]) 
write: 

a\ = t(Ai c cosh uj^x + A^ c cosh uj7x 



Af s sinhwf x + A ls sinha^ x) 



b\ = B^ c coshcj^x + B lc coshu^ x 
+ B^ s sinhw^x + B^ s smhuj7 x 

a,2 = — (Aic coshw2^ + Ais sinhw2^) 
k 

b-2 = B2C cosh 1^2% + -62s sinh L02 x 
a-3 = -j- exp(-w 3 x) 
63 = B 3 exp(-w 3 x) 



(8) 



with 



±\2 



p ± ak 

m 



P_ 



(9) 



(uj 3 ) 2 = k 2 + ^- and 5R(cj 3 ) > 



and where a 3 and 63 are the solutions in the region 3 for 
the non periodic problem only. In that case we applied the 
condition b3 — » when r — > 00. Furthermore, for sake of 
generality, we shall replace 03 by only in the numerical 
applications. 

Applying the appropriate symmetry conditions in x = 0, 
the following relations are found for the even solutions: 



B = 3?[b(x) exp(pt+ikz)] with b = (— ika, b, d x a), (5) 



-47 



.4 



Is 



B 



b: 



0. 



(10) 
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From ((6]) we have the additional relations: 
= B\ci A lc = —B lc 



with 



(11) 



for both problems non-periodic and periodic. For the even 
solutions of the periodic problem the boundary condition 
b 2 {R + 2e) = b2(R) leads to the additional relations 

s'mhuj2(R + e) = — A 2s coshu2(-R + e) 
B 2c smhuj2(R + e) = -B 2s cosh u> 2 (R + e). (12) 

For the odd solutions the same relations as [)10p and (jlip 
are found with inverted subscripts s and c. For the odd 
solutions of the periodic problem the boundary conditions 
b2{R + 2e) = — b2(R) leads to the additional relations 
given by (|12p with, again, inverted subscripts s and c. 



3.2 Boundary conditions 

The normal component of B, the tangential component 
of B//x and the z-component of the electric field E z — 
r/(V x B) z are continuous across each interface x = R 
and x — R + e. We can show that this set of relations is 
sufficient to describe all the boundary conditions of the 
problem. They write for 1=1, 2 (I = 1) corresponding to 
the non-periodic (periodic) case: 



-bi+i(x = xi) 



ai(x 


= x{) 




= Xl) 


—bi(x 


W 




= xi) 


— a[(x 


Vi 




= Xl) 



1 



Xl) 



J i+i 



(13) 



where the prime denotes the ^-derivative, x% = R and 
X2 = R + e. 

3.3 Resolution 

Replacing flgj) into |T3D and applying ([TD]). (fTTjl and (fH|) . 
we find a system to solve. The solution is non trivial only 
if the determinant of the system is equal to zero. This 
determinant writes 

1 + ¥+JL( X + + x -) +MSX + x - = (14 ) 

with X^ 1 = Rujf 1 tanh Rujf 1 for the even solutions and 
Xf~ = Ru>fcoth.Rujf for the odd solutions. 
For both even and odd solutions of the non-periodic prob- 
lem, M and S are defined by 

1+ m £ 3 tanhet^ s 1 + ™a tanhew 2 

M= Z " 2 : S = 



Ruj 2 ^+tanhew 2 



Ruj 2 ^+tanhew 2 



For the periodic problem, M and S are defined by 



(15) 



(M,S) = (m, s) C °^ 1 eUJ2 for even solutions, (16) 
R0J2 

(M,S) = (m, s) ta SUJ2 for odd solutions, (17) 
Ruj 2 



Mi' 



Va 
M3 : 



(72 



^3 ' 



(18) 



4 Results 

4.1 General remarks 

We can show that the growth rate p has no imaginary 
part as shown in Appendix A. An additional and simple 
argument is that, as the a-effect does not depend on z, 
there is no preferred way along the z-direction for a mag- 
netic wave to travel. Therefore, the marginal instability 
solution corresponds to p = 0. This leads to L02 — ^3 = k. 
The results are given in terms of the dimensionless quanti- 
ties to, n, s, r defined in (JTSJ) and of k = kR, p = pR 2 / 771, 
uji = Ru>i, e = e/R and R a = Rajr\\. The hat is dropped 
in the rest of the paper for sake of clarity. Then for a given 
set of the parameters m, n, s, r, k 7 e we seek solutions such 
that p — and R a is minimum. 

From © we see that replacing R a by -R a is the same as 
replacing k by — k and also the same as replacing ui^ by 
cj-f. As (|14p is symmetric in ui^ and , it is then suffi- 
cient to consider only positive values of k and R a . 
For the non-periodic problem we set 03 = (the region 3 
being insulating). This corresponds to the limit r — > 00 . 
It leaves (HI) unchanged but uj 3 = k and S = s tan ^ U2e . 



4.2 Influence of wall thickness 

In Fig. [3] the marginal curves \og w {R a ) versus log 10 (fc) 
are given for s = 10, to = n = 1 and different values of e. 
The marginal curves of the odd solutions are located be- 
tween the curves (d) and (e) and therefore always above 
the marginal curves of the even solutions located between 
the curves (a) and (c). Therefore only the even solutions 
are present at the onset of the dynamo action unless k is 
large. The same comments apply when varying s, to or 
n. Therefore in the subsequent subsections we shall focus 
only on the even solutions. 

The marginal curves of the non-periodic problem are lo- 
cated between the curves (c) and (b). They are always 
above those of the periodic problem located between the 
curves (a) and (b) . Then the periodic problem appears to 
be always more unstable than the non-periodic problem. 
However in the limit of large ke both problems have the 
same marginal curve (b). Indeed taking ke » 1 in (|15p 
and (|T6|) leads to (M,S) = (m,s)/k for both problems. 
Also taking ke >> 1 means that the fluid is embedded 
between walls of infinite thickness (compared to the ver- 
tical wave length of the field). Then no distinction can 
be found between both problems. From asymptotic esti- 
mates, we can show that 



lim 

ke — too.fe- 



R 



-1/2 



(19) 



which corresponds to the asymptotic left part of the curve 
(b). 



R. Avalos-Zuniga & F. Plunian: Influence of walls electromagnetic properties onto the onset of a stationary dynamo. 5 



+4 



logio(-Ra) 



o 



-2 









-5^H3r 






(b) 
















4 -2 





+2 



lo gioO) 



Fig. 3. The marginal curves log 10 (J? Q ) versus log 10 (fc) for m = 
n = 1, s = 10 and different values of e. The dashed curves 
between (a) and (b) correspond to the periodic even solutions. 
From (a) to (b): e=0 (curve (a)),e = 1, 10, 10 2 , 10 3 , 10 4 and 
e — > oo (curve (b)). The dotted-dashed curves between (b) 
and (c) correspond to the non-periodic even solutions. From 

(c) to (b): e—Q (curve (c)),e = 0.1, 1, 10, 10 2 and e — > oo (curve 
(b)). The dotted curves between (d) and (e) correspond to the 
odd solutions for both problems (periodic and non-periodic). 
They are obtained respectively in the limit of large e (curve 

(d) ) and for e = (curve (e)). 

The curve (a) is obtained for e = for the periodic prob- 
lem (even solutions). We can show that it is given by 

R a = k. (20) 

which has already been obtained in other periodic prob- 
lems (without walls) in which an anisotropic a-effect is 
the dynamo mechanism (see e.g. |16j|18j). 
The curve (c) is obtained for e = for the non-periodic 
problem (even solutions). From asymptotic estimate we 
can show that for k < < 1 it is given by 

R a = y/3n/mk. (21) 

In the limit of small or large values of both k and ke, 
asymptotic expansions in k of R a have been calculated. 
A synthesis of these expansions is given in table [TJ These 
expansions fit very well to the slopes of the curves of Fig. 

m 

In the limit of small ke and large e these expansions also 
give R a = 0(e -1 / 2 ) for the non-periodic and R a = O(e) 
for the periodic even solutions. Then for increasing e, the 
non-periodic problem gets more unstable whereas the pe- 
riodic problem gets more stable. To get a qualitative ex- 
planation to this striking difference between both prob- 
lems, we sketch in Fig. [4] the dimensionless dissipation 
rate J (ratio of the Joule dissipation to the magnetic en- 
ergy, see details in Appendix B) versus e. We also sketch 
on the same figure the dimensionless rate of the work of 
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lo Sio( e ) 

Fig. 4. Dissipation rate J and work of the Lorentz forces rate 
5 versus log 10 (e) for m = n = 1, s = 10 and k = 0.01 for both 
problems periodic (P) and non-periodic (NP). 

the Lorentz forces S. Then R a is defined by R a = J jS. 
What is remarkable is that the decrease (increase) of R a 
for the non-periodic (periodic) problem is not only due 
to a decrease (increase) of J but also to an increase (de- 
crease) of S. We also find that the dissipation in the wall 
in both cases is always negligible compared to the dissi- 
pation in the fluid. This shows that changing e leads to 
a pure geometrical effect on the field and current lines, 
a consequence of it being the change of the dissipation 
in the fluid and work of the Lorentz forces. As an illus- 
tration we sketch in Fig[5]the isolincs of $t[b(x) exp(ikz)] 
for different values of e. These isolines correspond to the 
current density lines in the (x, z)-plane and also to the iso- 
lines of the y-component of the magnetic field. We see that 
increasing the thickness (from left to right) for both prob- 
lems periodic (top) and non-periodic (bottom) changes 
indeed the geometry of these isolines. Increasing e leads 
to a bending (flattening) of the current lines for the peri- 
odic (non-periodic) problem, enhancing (decreasing) Joule 
dissipation, consistently with Fig[U 

4.3 Influence of wall conductivity and permeability 

The expansions of table [T] for small values of both k and 
ke also give the dependency of R a versus m and s. For 
large (resp. small) value of m and s, R a = O^ms)^ 1 ^ 2 ) 
for the non-periodic (resp. periodic) even solutions. This 
means that the dynamo onset is easier to reach with walls 
of higher permeability and/or conductivity. This is true 
also outside the previous asymptotic limits as depicted in 
Fig. H Fig. E| and Fig. El 

In Fig.[6]the marginal curves R a versus log 10 (s) are given 
for fc = l,m = n= l and different values of e. The dashed 
(dotted-dashed) curves correspond to the periodic (non- 
periodic) even solutions. Again, both problems (periodic 
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/7r\2 /o m-f s-f 2\/2ms 1 I r\( 1 ~\ 
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T + *0(1) 


fc + ^ - v2'-r + + s 2 /jr s h) + o(^) 



Table 1. Asymptotic expansion in k of R a for both problems non-periodic and periodic. The parameter c (-| < c < 7r) is solution 
of the equation ^ (tan c + tanh c) = — c for the non-periodic problem and + msec cot c)( + msec coth c) = ( lr ^) 2 for 

the periodic problem. 
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Fig. 5. Isolines of $l[b(x) exp(ifcz)] in the (a;, z)-plane (0 < 
x < 1) for m = n = 1, s = 10, fc = 0.01, for both periodic 
(0 < kz < 7r/100) and non-periodic (0 < kz < n) problems 
and different values of e. The isolines of the periodic problem 
are given for (a) e = 10 _1 , (b) e = 1, (c) e = 10 and (d) 
e = 10 2 . The isolines of the non-periodic problem are given for 
(e) e = 0, (f) e = 10" 2 , (g) e = 10 _1 and (h) e = 1. 



and non-periodic) have a common solution corresponding 
to the full curve of Fig. [6] in the limit of large e. The 
monotonic decrease of R a versus s is in agreement with 
the stationary solutions of [3] for the outer wall problem 
and of [12] for the inner wall problem. 
In Fig.[7]the marginal curves R a versus log 10 (m) are given 
for fc = l,s = n = l and different values of e. The dashed 
(dotted-dashed) curves correspond to the periodic (non- 
periodic) even solutions. Again, both problems (periodic 
and non-periodic) have a common solution corresponding 
to the full curve of Fig. [7] in the limit of large e. This 
curve coincides with the full curve of Fig. [6] Indeed for a 
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Fig. 6. The marginal curves R a versus log 10 (s) for fc = m — 
n = 1 and different values of e. The dashed (dotted-dashed) 
curves correspond to the periodic (non-periodic) even solu- 
tions. The labels indicate log 10 (e). The full line is common 
to both problems in the limit of large e. 



given k and in the limit of large e we find from (|15[) and 
(fT6| that (M, S) = [m,s)/k for both problems. As dH]) 
is symmetric in M and S 1 , the full curves of Fig. [S] and 
Fig. [7] coincide. For the same reason the dashed curves 
corresponding to the periodic problem of Fig. [6] and Fig. 
[7J coincide. Again the monotonic decrease of R a versus m 
is in agreement with the stationary solutions of [3] 
In Fig. [H]the isolines of R a are sketched in the (log 10 ( s ); 
log 10 (m))-plane for k — 1, e = 0.1 and n — ra. The iso- 
lines of the periodic even solutions (dashed curves) are 
symmetric to the straight line m = s from the same sym- 
metry arguments as above. This is not true for the non- 
periodic even solutions (full lines). In that case the asym- 
metry comes from the fact that 173 and /i3 do not play a 
symmetric role as 03 = and /i3 7^ 0. In table [2] we give 
the values of R a for asymptotic values of (s,m). We note 
that they do not depend on the type of the problem either 
periodic or non periodic. 

At that point we can make two conclusions. First in the 
non periodic case when varying the electromagnetic prop- 
erties of the outer wall, the roles of the conductivity and 
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Fig. 7. The marginal curves R a versus log 10 (m) for k = s = 1, 
n = m and different values of e. The dashed (dotted-dashed) 
curves correspond to the periodic (non-periodic) even solu- 
tions. The labels indicate log 10 (e). The full line is common 
to both problems in the limit of large e. 
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Fig. 8. Isolines of R a in the (log 10 (s), log 10 (m))-plane for 
k = 1, e = 0.1 and n — m. The dashed (full) curves correspond 
to the periodic (non-periodic) even solutions. 





s « 1 


s » 1 


771 << 1 


3.5 


2.04 


m » 1 


2.03 


1 



Table 2. Values of R a for asymptotic values of (s, m) and 
k — 1, e — 0.1, n = m and both problems periodic and non 
periodic. 



permeability are different. In other words, a given jump 
of magnetic diffusivity n = X/afj, between the fluid and 
the wall can lead to different results depending on the re- 
spective jumps of the conductivity and permeability. Sec- 
ond, for the periodic problem and e = 10 _1 , we see from 
Fig. [7] that when increasing m from 1 to 10 3 the reduc- 
tion of the threshold is small (4.5%). Now we see from 
Fig. [5] that such a reduction at m = 10 3 can be compen- 
sated by the decrease of s from 1 to 0.89. In other words 
the threshold is the same (<~ 1.1) for (s,m) = (1, 1) and 
(s, m) = (0.89, 10 3 ). Then the reduction of the threshold 
when using ferromagnetic steel inner walls (with a relative 
conductivity of order 1) instead of non ferromagnetic ma- 
terials is negligible. Therefore it does not seem likely that 
the dynamo action in a FBR can be favored by the use of 
ferromagnetic assemblies at least in the framework of this 
simplified study. This is in disagreement with more elab- 
orate models [17] which show that the jump of m from 1 
to 10 3 reduces the threshold significantly enough to start 
the dynamo action inside the core of a FBR. A reason 
of discrepancy may come from the fact that in [17] s is 
kept constant. Also it may come from the idealized ge- 
ometry of the flow inside each assembly in [17\ which is 
helical instead of being replaced by some a-effect as here. 
Now let us consider the influence of a ferromagnetic belt 
surrounding the core of a FBR. From Fig. [7] we see that, 
for the non periodic problem, from m = 1 to m = 10 3 the 
threshold reduction rate is significant (larger than 20% for 
e = 10 _1 and e = 1). If now we assume that the ferromag- 
netic material is less conducting than the fluid then from 
Fig. [8] jumping from (s,m) = (1, 1) to s < 1 and m = 10 3 
leads to a minimum threshold reduction of 10%. Finally 
it is likely that an outer ferromagnetic belt can favour the 
dynamo action in the core of a FBR. 



4.4 Influence of the fluid permeability 

The influence of the fluid permeability can also be given 
by the expansions of table [T] for small values of both k 
and fee. We see indeed that R a = O^n/m) 1 / 2 ) for the 
non-periodic even solutions and that R a = (^((l/m) 1 / 2 ) 
in the limit of large 1/m for the periodic even solutions. 
Therefore for both problems a fluid with large permeabil- 
ity leads to higher R a . This is also confirmed by Fig[9]in 
which the marginal curves R a versus n/m are given for 
fc = l,s = n = l and different values of e. The dashed 
(dotted-dashed) curves correspond to the periodic (non- 
periodic) even solutions. Though R a increases versus n/m 
we must note however that R a jo\\i\ decreases versus n/m 
(see also [3]). This means that, keeping R a constant and 
using a fluid of high permeability allows the experimenter 
to reduce significantly the flow intensity or the size of the 
device. Assuming as in [21j that the power to drive an 
experiment is dissipated by turbulence, leads to a power 
proportional to U 3 . Therefore using a fluid of high perme- 
ability would imply a significant reduction of the driving 
power. This has also motivated some experimental studies 
[22] [23]. 
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Fig. 9. The marginal curves R a versus log 10 (n/m) for k = 1, 
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or 

nxB=jn = 0. (23) 

The first boundary condition (j2"2"|) corresponds to region 
2 being a perfect conductor whereas the second boundary 
condition (|23|) corresponds to region 2 having an infinite 
permeability. A priori we would expect to recover these 
two limits with our model taking the limit s » I or 
m » 1. Furthermore these limits should not depend on 
e nor on the type of problem periodic or non periodic 
which is considered. As can be seen from Fig. [6] and Fig. 
[7] this is not true. In fact it is known [Mj that considering 
an ordinary body of electric conductivity a and magnetic 
permeability //, the supraconductivity limit corresponds 
not only to a — > oo but in addition to /J, — ► . Considering 
the double limit s » 1 and m << 1 in our model we 
indeed recover the boundary condition (|22|) . A body of 
high permeability is also known to be a poor conductor 
and corresponds not only to \i — » oo but in addition to 
a — > . Considering the double limit m > > 1 and s < < 1 
in our model we indeed recover the boundary condition 
(|23|) . For both double limits the dynamo threshold is found 
to be R a w 2 as given in tableau- 





e/ 


»/ 


m = n /" 


1/m/ 


Non per. 


5 


/ 7.7 


/ 1.17 


\ 6.22 


/ 6.22 






\ 6.76 


\ 1.65 


\ 11 


/ 11 


Ra = 




\ 52 


\ 1.93 


\ 1.77 


^ 1.77 


Periodic 


5 


\2 


<" 1.14 


\ 1.3 


/ 1.3 






/ 17 


\ 1.77 


\ 2.52 


/ 2.52 


Ra = 


J/S 


\ 34 


\2 


\ 1.94 


' 1.94 



Table 3. Global information on the evolution of S and J7 
versus e (for fc = 0.01) and versus s, m — n and 1/m (for 
fc = 1). The arrows indicate if the quantity grows or decays. 
The numbers indicate the factor of growth or decay. 



4.5 Geometrical effects 

Following the same idea as in section 14.21 we calculate the 
dissipation and work of the Lorentz force rates for each 
previous case (changing the conductivity or permeability 
of the wall or the fluid) . We find again that the change of 
J is always accompanied by a change of S and that the 
dissipation in the wall is always negligible compared to 
the dissipation in the fluid. In tabled we give some global 
informations on the behavior of J and S for the different 
previous cases. For the same reasons as in section T4. 2 1 we 
conclude that adding outer or inner walls with different 
conductivity or permeability has a pure geometrical effect 
on the field and current lines. 

4.6 Infinite conductivity or permeability 

Two usual ways to simplify the dynamo problem are to 
solve the induction equation in region 1 only, with one of 
the following boundary conditions at x = 1: 

B-n = nxE = 0, (22) 



5 Conclusions 

For a stationary dynamo with either insulating or periodic 
boundary conditions, we showed that additional outer or 
inner walls change the geometry of the field and current 
lines in the fluid. This geometrical effect has an effect on 
the dynamo instability threshold. For given inner or outer 
walls, increasing their conductivity or permeability helps 
for dynamo action. In the other hand increasing the thick- 
ness of the inner wall plays against the dynamo action 
(contrary to the outer wall for which increasing the thick- 
ness helps for dynamo action) . These conclusions are con- 
sistent with those obtained in other geometries for outer 
3J [8] [9] as well as inner [12] walls and we believe that they 
are generic in the sense that they do not depend on the 
generation process (a 2 or else) as far as the solution stays 
stationary. The detailed mechanism of the geometrical ef- 
fect is however non trivial as the dissipation in the fluid is 
changed as well as the work of the Lorentz forces. In any 
case the dissipation in the walls is always negligible. 
The usual boundary conditions used to describe a per- 
fectly conducting or high permeability outer wall are re- 
covered with a double limit on s and m, stressing that a 
supraconductor outer wall would correspond to s — > oo 
and ttl — ► and a high permeability outer wall to m — > oo 
and s — > 0. 

As mentioned in section 12.11 an additional mean flow U 
like in the core of a fast breeder reactor might lead to os- 
cillatory dynamo solutions. Then in this case as shown in 
[3] for the Riga dynamo experiment, some additional eddy 
currents in the walls may lead to significant dissipation in 
the walls reminiscent to a skin effect. Then the previous 

1 We believe that the fact that R a is the same for both double 
limits is coincidental and probably related to our model. 
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behavior of R a versus e, s or m would not be monotonic 
anymore. For the specific case of the core of a FBR, the 
conclusions of section 14.31 should then be revised with re- 
spect to this additional mean flow U. It is not clear a priori 
what would be the main effect of U, either increasing or 
decreasing the dynamo threshold. 

R. Avalos-Zufiiga acknowledges the Mexican CONACYT 
for financial support. 

Appendix A. Proof D that p has no imaginary 
part 

Taking (J7JI and the complex conjugate of ©, we easily 
find 

rjiaig// - (i]ik 2 + p)\ai\ 2 - a^ai = (A.l) 
mb'i% - (riik 2 +p)\k\ 2 - ai k 2 b iai = (A.2) 

where an underlined quantity means the complex conju- 
gate of this quantity. Integrating by part and using the 
boundary conditions (|13[) we can show that 

m 



aigf/dx = az-iaj_i(£;-i) +a/aj(x;) 



b'fadx = - 



\a\\ 2 dx 

xi-l 

Vi-UM u 



(A.3) 



\b[\ 2 dx 



(A.4) 



where xq — ao = bo = 0. For the non-periodic case 
X3 — +oo and a^(x — +oo) = b^(x = +oo) = 0. For 
the periodic case the periodicity implies that a,2g/ 2 (x = 
R + e) = b' 2 b 2 (x = R + e) = 0. The previous relations 
do not depend on the parity of di nor b±. Combining the 
integral of relations (|A. 1[) and (|A.2|) . we find: 



° = J2- T (-mk'latf + mlbtf-mk 4 ^ 2 



ai f (A.5) 



+ m k 2 \bi\ 2 -^-k 2 p\ ai \ 2 +p\b l \ 2 )dx. 

m - 



Then it is straightforward to show that = 0. 



Appendix B. Expression of the dissipation rate 

Multiplying (|4]) by B//i we obtain the following equation 
(where underlying means complex conjugate): 

|(^) = f(^)-«(b-(z.b)z).j-^ (B.l) 
at 2/i ox [i -a 



2 A similar proof of the principle of exchange of stabilities 
has been given for an a 2 -dynamo with a constant a within an 
electrically conducting sphere surrounded by an insulator [25] . 



with b = (—ika, b, a'), j = j^(~ikb, —Aa, b'), e = j/cr + 
a(b — (z.b)z). Applying the boundary conditions (| 13[) we 
can show that J Q F -g-{=^)dx = where r = R + e (= 
+oo) for the periodic (non-periodic) problem. At the dy- 
namo onset J^(|b| 2 /2/i) = and the Joule dissipation 
|j| 2 /cr is equal to the alpha-dynamo source — a(b — (z.b)z)- 
j. In that case (at the onset) the expressions of the dimen- 
sionless dissipation J and work of the Lorentz forces S 
rates are given by 



J 



ti + "l{{k 2 + ^)\b\ 2 + \V?)dx 

Jo ±(k 2 \a\ 2 + \b\2 + \ai\2)dx ' 
Jo j(ba" -2k 2 ab)dx 

f o ^(fc2| a |2 + |5|2 + H 2 )cK ' 



(B.2) 



(B.3) 



where x = x/R, f\ = rj/r]i, jl = fi/fii, k = Rk and b — Kb. 
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